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Context. Dwarf galaxies can arise from self- gravitating structures emerging from tidal tails. What fraction of the known dwarf galaxies 
in the Local Universe can have this origin is still a matter of debate. 

Aims. In our elfort to understand the origin and evolution of tidal dwarf galaxies and their correspondence with local objects, the first 
step is to understand how these galaxies (which are supposed to have a limited amount of dark matter) react to the feedback of the 
ongoing star formation. 

Methods. We make use of 2-D chemodynamical calculations in order to study the early evolution of isolated, dark matter-free dwarf 
galaxies. We present models in which feedback parameters are varied. We also compare the results with dark matter-dominated dwarf 
galaxy models. 

Results. All the considered models show that the star formation proceeds for more than 300 Myr, therefore dwarf galaxies without 
large dark matter halos are not necessarily quickly destroyed. The chemical evolution of these objects is consistent with the main 
chemical properties of the dSphs of the Local Group. Models with large dark matter halos show results consistent with models free 
of dark matter, indicating that the distribution of gas is more important than the depth of the potential well in determining the global 
behaviour of dSph-sized dwarf galaxies. 
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1. Introduction 

Galactic collisions and close encounters are thought to be com- 
mon events in the Local Universe. It has been clear for some 
time that the structure of a galaxy can be affected by tidal 
. !^ interactions with close neighbours (Toomre & Toomre |19721 l. 
Even for what concerns the Milky Way, signs of the interaction 
$_i with the Magellanic Clouds have been analysed in detail (e.g. 
5^ Mastropietro et al. 2005; Ruzicka et al.'TQQT). 

There is growing evidence that tidal arms, produced by the 
interaction of gas-rich galaxies, can trigger the formation of self- 
gravitating structures evolving into objects similar to the local 
dwarf galaxies (Due & Mirabel 1998; Walter et al. 2006 ). These 
objects are therefore called tidal dwarf galaxies (TDGs) and 
their production rate must have been higher in the early universe 
than today, due to the higher probability of tidal interactions and 
the larger gas content of young galaxies. From a stellar dynam- 
ical point of view, it is shown that these star complexes typi- 
cally survive the first few hundred Myr (Kroupa [l9981 l and can 
evolve into objects similar to the present-day dwarf spheroidal 
(dSph) galaxies (Kroupa [T997l Metz & Kroupa llOOTl l. Recently, 
detailed numerical simulations of the evolution of tidal tails have 
been presented (Bournaud & Due 120061 Wetzstein et al. 2007 ) 
showing that the formation of TDGs in these sites is favoured 
by the presence of large and extended reservoirs of gas. These 
studies also show that TDGs must have a very reduced dark 
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matter (DM) content, confirming previous results (e.g. Barnes 
& Hernquist 1992 ). Okazaki & Taniguchi ( 2000 ) even point out 
that the majority, if not all the dwarf galaxies could be TDGs 
if only a few long-Uved TDGs emanate from each late-type en- 
counter. 

It is therefore important to understand the early evolution of 
TDGs, in order to clarify how the potentially disruptive feed- 
back from the ongoing star formation (SF) acts on a DM-poor 
structure, characterised by a very reduced potential well, and to 
compare them with DM-rich models. Finally, the study of the 
chemical evolution of dwarf galaxies and comparison with ob- 
served ones offers a wealth of information about the past history 
and the origin of these objects. 

In this paper we study for the first time the 2-D chemody- 
namical evolution of DM-free model galaxies, in order to under- 
stand under which conditions these objects can sustain the en- 
ergy released by the dying stars without experiencing a complete 
blow-away, as well as what the final chemical patterns are. This 
work is an extension of the 1-D chemodynamical simulations 
of DM-free dwarf galaxies of Hensler et al. (2004). In Sect. |2] 
we describe the adopted model and in Sect. [3] we summarise the 
main results of our investigations. Finally, in Sect. |4] we draw 
some conclusions. 



2. The model 

We performed 2-D simulations of the dynamical and chemical 
evolution of TDGs in cylindrical coordinates. The hydro solver 
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and the implementation of the routines for the chemical evolu- 
tion have been described in Recchi et al. (' 200 U [2004 1) and refer- 
ences therein. Several new implementations have been added in 
this version and will be comprehensively described in Recchi et 
al. (in preparation). The main code extensions are: 

- self-gravity like in Rieschick & Hensler (I2003I I. which solves 
the Poisson equation by means of the so-called "Alternate 
Direction Implicit" method. 

- Metal-dependent stellar wind luminosities taken from the 
Starburst99 (Leitherer et al. 1999) software package. 

- A star formation rate proportional, within each grid cell, to 
the amount of gas (through a constant ssf) with an upper 
temperature threshold (Tthr) and valid only for convergent 
gas flows (i.e. when V ■ v < 0). ssf denotes the fraction of 
gas converted to stars within the next 5 Myrs, once the SF 
criteria are fulfilled. 

The calculation strategy is that: the SF criteria identify the 
grid points that, within an interval of 5 Myr, produce stars. The 
star masses are distributed according to a Salpeter IMF. All the 
relevant variables concerning this SF generation are stored in ar- 
rays that are used to calculate the feedback from SNela, SNell, 
and intermediate-mass stars. In this way, the chemical composi- 
tion of stars is always recorded and metal-dependent yields and 
energy production can be taken into consideration properly. In 
spite of some tricks to save memory, the calculations are very 
computationally demanding (several weeks to months for each 
run), therefore only a limited number of runs can be performed. 
For our convenience, the dynamics of the stars is neglected in 
this set of models, therefore the stars remain where they are born. 
Due to their short lifetimes, this is not critical at all for the mas- 
sive stars, which are the main sources of stellar energy feedback. 

Although these simulations are not aimed at reproducing 
specific objects, the Local Group dwarf galaxy Fornax gives us 
a good reference for calibrating the initial setup. Indeed, it has 
been suggested that the dSphs near to the Milky Way may be 
TDGs because their spatial distribution is incompatible with a 
cosmological origin (Kroupa et al. 120051 ). Although Local Group 
dSphs are commonly believed to be DM-dominated, there are 
arguments suggesting that they are not necessarily DM-rich. 
Twisted and squashed isophotes are common in Local Group 
dSphs and are not compatible with the presence of a massive 
DM halo. Most of the kinematical properties of these objects 
can also be explained without DM (Metz & Kroupa [20071 ) . An 
extensive discussion about this point will be presented in Recchi 
et al. (in preparation). The initial gas distribution is spherically 
symmetric and follows a King profile with a core radius of 500 
pc. The total gas mass in the grid (a sphere of ~ 8 kpc of radius) 
is calculated according to Mg - Mpomax / ssf, Mpomax being the 
mass of Fornax (assumed to be 6.8 ■ 10^ Mq, Mateo |T998l ). 

We consider models in which the SF efficiency ssf and the 
temperature threshold Tthr are varied. Two models with a large 
rigid DM halo are considered, too, in order to outline the dif- 
ferences with DM-free models. The central DM density is cal- 
culated according to the observed velocity dispersion of Fornax 
(10.5 km s"', Mateo 1998) and the distribution follows a quasi- 
isothermal profile with core radii of 1 kpc and 200 pc, respec- 
tively. These models contain 1.1 ■ 10^ and 3.7 ■ 10^ Mq of DM 
within 1 kpc, respectively. We also consider a model in which 
the feedback from stellar winds is suppressed. Model parame- 
ters are summarised in Table (Tj A wider range of values of ssf 
and Tthr have been considered, as well as models with better spa- 
tial resolution and models with an initial anisotropic distribution 



of gas. For the sake of simplicity, we focus in this letter only 
on the 6 basic models, mentioning briefly the behaviour of the 
other models. Only the first 300 Myr of the evolution have been 
calculated. This is due to the large computational time required 
by the simulations, but also to the fact that for timescales longer 
than this, interaction effects with surrounding galaxies are sup- 
posed to play a major role in the evolution of these small objects. 
The effect of tidal interactions on the evolution of these objects 
is currently under investigation. 

Table 1. Model parameters 



Model 


DM? (core radius) 


£SF 


T,h, (K) 


Stellar winds? 


1 


No 


0.1 


10* 


Yes 


Idc 


Yes (200 pc) 


0.1 


lo-* 


Yes 


Ide 


Yes (1 kpc) 


0.1 


lo-* 


Yes 


2 


No 


0.1 


5 • 10' 


Yes 


3 


No 


0.1 


lO'* 


No 


4 


No 


0.05 


lO'* 


Yes 



3. Results 

3.1. A typical evolutionary sequence 

We show here an example of the chemodynamical evolution 
of an isolated, DM-free dwarf galaxy. This model (designated 
as Model 1, see Table [TJ is characterised by esF = 0.1 and 
Tthr = 10"^ K. This temperature is much higher than the typical 
temperatures of the cores of molecular clouds, where SF takes 
place. Indeed the size of each computational cell is also much 
larger than a molecular cloud core. (The size of the first cell is 
10 pc, then it proceeds radially outwards with a size ratio be- 
tween adjacent zones of 1.028.) Since this temperature has to be 
considered as the average temperature of a region of the galaxy 
encompassing a star-forming molecular cloud, Tthr is necessarily 
a free parameter in our models. We discuss its effect in Sect. 13.21 
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Fig. 1. Density contours of the gas for Model 1 at 9 evolutionary 
times (labelled, in Myr, at the top right corner of each panel). The 
density scale (in g cm"^) is on the righthand strip. 



The evolution of Model 1 during the first 150 Myr is shown 
in Fig. [T] During the first few 10^ yr, the central density is not 
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high enough to trigger SF; therefore, the system is only charac- 
terised by a slow infall motion, due to the self-gravity. After ~ 30 
Myr the gas has been condensed enough to form the first stars. A 
network of cavities with increasing size and complexity is cre- 
ated due to the patchy distribution of the SF sites. It maintains, 
however, a spherical symmetry and creates a gradually expand- 
ing supershell. The superbubble grows slowly in size, reaching 
~ 600 pc in 150 Myr Only at t ~ 200 Myr does the size of the 
superbubble reach 1 kpc, and at this point the SF rate (SFR) be- 
gins to decrease, due to the reduction of available gas. At the end 
of the simulation, it is approximately one fourth of its peak value 
(see Fig. |2|l. The evolution of gas mass with time in the central 
kpc of Model 1 (Fig.|2l third panel) is very similar to the evolu- 
tion of the SFR with time. This is not a surprise, since we have 
chosen a linear dependence between SF and the gas mass. The 
increase in the gas mass in the first ~ 150 Myr is due to the fact 
that the density distribution extends up to several kpc and the 
self-gravity attracts a fraction of the gas located in the external 
regions. The global metallicity of the stellar component (Fig. |2l 
lower panel) oscillates around a mean value of ~ 0. 1 Z©, consis- 
tent with what is observed in Fornax and in the largest dSphs of 
the Local Group (Mateo IT998I ). For this model we checked the 
timescale required to blow-away the gas from the central kpc, 
and this corresponds to ~ 400 Myr We reran this model with 
higher spatial resolution, i.e. a resolution of 5 pc. This model 
shows a very similar behaviour to Model 1 . 
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t (Myr) 

Fig. 2. Evolution with time of the SFR (in M© / yr, upper panel), 
stellar mass (in 10^ Mq, second panel), gas mass within 1 kpc (in 
10** Mq, third panel), and mean metallicity of the stellar compo- 
nent (in Zq, lower panel) for Model 1 (reference model, very 
thick lines). Model 2 (reduced Tthr, thick lines). Model 3 (no 
stellar winds, thin lines), and Model 4 (reduced ssp, very thin 
lines). 



With our model, we are also able to check the build-up of the 
[O/Fe] vs. [Fe/H] relation in the stellar population. This is shown 
in Fig. [3] The evolution in the first few 10^ yr is dominated by 
SNell ejecta and this produces a plateau of [O/Fe]. At [Fe/H] ~ - 
2, the Fe produced by the SNela (which explode with some delay 
compared to SNell, Matteucci & GreggiO [19861 l starts reducing 
the [O/Fe]. For Local Group dSphs, a knee at [Fe/H] ~ -2 (much 
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Fig. 3. [O/Fe] vs. [Fe/H] in stars. Symbols as in Fig.|2] 



earlier than the analogous knee observed in Milky Way stars) is 
a common feature (see e.g. Tolstoy et al. i2003. Lanfranchi & 
Matteucci |2004] l. 

3.2. A parametric study 

With its lower temperature threshold for the onset of SF, Model 2 
is supposed to be characterised by a reduced impact of the feed- 
back and therefore a milder SFR due to the fact that a smaller 
amount of gas fulfils the SF condition. Indeed, its behaviour 
(Figs. |2] and O is similar to the reference model, showing that 
SF self -regulation (Koppen et al.'1995) is at work, similar to that 
in 1-D chemodynamical modelling of dwarf galaxies (Hensler et 
al. l2004l l. 

In Model 3 we have suppressed the input of energy com- 
ing from stellar winds, thereby reducing the feedback. There is 
a significant effect on the global thermal energy budget, result- 
ing in a slightly larger amount of gas inside the galaxy at the 
end of the simulation and, consequently, in a larger SFR. This 
model is more stable against the feedback of the dying stars, so 
a longer duration of SF is expected. At the end of the simulation, 
~ 10^ Mo of stars are already present inside the galaxy (Fig.|2]i, a 
higher value than the estimated stellar mass of Fornax. We have 
to take into consideration, however, that more than 60 % of these 
stars will die within 10 Gyr. Moreover, tidal stripping is likely to 
produce a substantial reduction of the stellar component (Kroupa 
[1997 ). 

Model 4, with a reduced ssf, has a higher initial gas mass 
(see Sect. 2) and therefore starts forming stars earlier and at 
a higher rate. Its attained metallicity is of the order of a few 
hundredth of Z©, considerably lower than the observed value in 
Fornax. This is mostly due to how the newly produced metals are 
diluted in a larger fraction of pristine gas but also to the fact that 
the reduced injection energy rate per unit mass into the system 
reduces the turbulence, therefore decreasing the diffusivity of the 
medium and increasing the mixing timescale. This value of the 
global metallicity, although inconsistent with the chemical com- 
position of Fornax, is similar to the metallicity of smaller dSphs, 
supporting the idea that the SF efficiency increases with the mass 
of the object (e.g. Lanfranchi & Matteucci 2004). The evolution 
of [O/Fe] vs. [Fe/H] for this model shows a large loop (Fig. [3), 
caused by the continuous increase in total gas mass through gas 
inflow between ~ 100 and ~ 200 Myr that also produces a de- 
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crease in the global metallicity (see Fig. |2]i. Therefore, at vari- 
ance with the other models, the [Fe/H] is not a proxy of the time. 

Two DM-rich models (models Idc and Ide) are also con- 
sidered. Their setup is identical to model 1, but their gravity is 
dominated by a large static DM halo. The comparison of the SFR 
in these models is shown in Fig.|4] The differences between the 
different models are very small. Also the other properties of the 
models are not significantly affected by the presence of a large 
DM halo. This shows that the global behaviour of the gas in the 
central part of the galaxy is mostly determined by the total mass 
and distribution of gas itself rather than by the depth of the po- 
tential well. However, it is foreseeable that the late evolution of 
the gas is affected by the presence of a DM halo extending up to 
very large distances. 




t (Myr) 

Fig. 4. Evolution with time of the SFR (in M© / yr), for Model 
1 (reference model, very thick line). Model Ide (extended DM 
halo, thick Une), and Model Idc (concentrated DM halo, thin 
line). 



We also ran a simulation in which esF=0.2. In this case, 
the expanding cavity created by overlapping SNe is energetic 
enough to produce a blow-away and unbind most of the gas in 
the galaxy in ~ 100 Myr This determines the duration of the SF 
period, which is very short and inconsistent with the typical SF 
histories of Local Group dSphs (Grebel [l997l) . 

Finally, when an asphericity in the initial distribution of gas 
is introduced, the gas tends to expand preferentially along one 
axis, namely the one in which the pressure gradient is steeper. 
Interestingly, the development of a polar galactic wind does not 
suppress the SF process, since some gas is compressed towards 
the walls of the galactic chimney where the condition for the 
onset of SF can be fulfilled (see Recchi et al. l2007l l. 

4. Discussion and conclusions 

By means of a detailed 2-D chemodynamical code, we have sim- 
ulated the evolution of isolated dwarf spheroidal galaxies with- 
out DM in order to study the early evolution of TDGs and com- 
pared the results with DM-dominated models. 

We have seen that reasonable assumptions about the SF effi- 
ciency and the temperature threshold above which the SF is sup- 
pressed lead to a SF lasting at least 400 Myr. This interval of time 
broadens further if we consider very small SF efficiencies (which 



is possible in small galaxies). Provided that it keeps its spherical 
symmetry, a small galaxy is therefore relatively stable against 
the feedback of the ongoing star formation, even in absence of a 
large DM halo. This issue has already been shown by Hensler et 
al. (2004). However, we have seen that a DM-dominated model 
shows a behaviour similar to the DM-free ones, indicating that 
the gas dynamics of the galaxy, at least in its central part, is dom- 
inated by the gas density and distribution more than by the depth 
of the potential well. Eventually, after several hundred Myr, a 
strong galactic wind may arise, able to eject most of the gas. 
However, this gas does not necessarily leave the parent galaxy. 
It may indeed later be reaccreted by the gravitational pull of the 
galaxy if tidal interactions do not disperse it. Recchi & Hensler 
(120061) have done experiments about the timescale required to re- 
fill with cold gas the centre of a galaxy after a supernova-driven 
galactic wind has occurred. This timescale is of the order of a 
few hundred Myr, after which a new episode of SF can occur 
This might lead to the bouncing back and forth of the gas, due 
to the competing effects of the gravitational pull and the feed- 
back of the ongoing star formation, resulting in a sort of self- 
regulation, until the moment in which a close encounter with a 
more massive galaxy definitely moves the gas away, shutting off 
the cycle. This mechanism is under investigation at the present 
time. 

A more detailed description of the code and a more extensive 
analysis of the parameter space will be provided in a paper in 
preparation. Further work and modelling is necessary to sharpen 
our understanding of TDGs. In particular we are including stellar 
dynamics and a time-dependent tidal field due to the interacting 
system. Also the IMF is going to be subject of detailed analysis. 
In particular, we will implement the integrated-galactic initial 
mass function (Weidner & Kroupa [2005l l. in order to consider the 
variations of galaxy-wide IMF as a function of the SFR. Finally, 
we have started 3-D simulations of TDGs (Marcolini et al., in 
preparation). 
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